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The critical Ising model in two dimensions with a defect line is analyzed to deliver the first 
exact solution with twisted boundary conditions. We derive exact expressions for the eigenvalues 
of the transfer matrix and obtain analytically the partition function and the asymptotic expansions 
of the free energy and inverse correlation lengths for an infinitely long cylinder of circumference 
Lx- We find that finite-size corrections to scaling are of the form aklL‘^~^ for the free energy / 
and hk{p )!and Ck{p )!for inverse correlation lengths and respectively, with 

integer values of /c. By exact evaluation we find that the amplitude ratios bk{p)/ak and Ck{p)/ak 
are universal and verify this universal behavior using a perturbative conformal approach. 
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The study of boundary conditions (EC’s) in confor¬ 
mal field theories has attracted much attention in recent 
decades. This is because of its relevance in string and 
brane theory on the one hand, and in various problems 
of statistical mechanics and condensed matter physics on 
the other. A classification of conformal boundary condi¬ 
tions on the torus has been given in [4]-[6]. Computation 
of the partition function involves identifying the states at 
the two ends of a cylinder through the trace operation. 
One may insert a defect line (or seam) into such a system, 
along a non-contractible circle at the end of the cylinder 
before closing it into a torus, the effect of which is to 
twist the boundary conditions. It is important to under¬ 
stand the effects of such a twist, and to do this it is valu¬ 
able to study model systems, especially those which have 
exact results. The Ising model is the most prominent ex¬ 
ample and one of the best studied models of statistical 
mechanics. Since Onsager obtained the exact solution of 
the two-dimensional Ising model with cylindrical EC’s in 
1944 [7], there have been continuous attempts to treat 
different two-dimensional topologies |8HI3]. 

At the critical point Tc the asymptotic finite-size scal¬ 
ing behavior of the critical free energy /tv and the inverse 
correlation lengths of an infinitely long 2D cylinder 
of finite circumference N has the form [14] 

lim [N^ifN - /oo)] = A, lim = Dn 

where /oo is the bulk free energy and A and Dn are the 
universal constants which may depend on the EC’s. The 
values of A and Dn are known to be related to the confor¬ 
mal anomaly c, the conformal weight A, and the scaling 
dimension of the n-th scaling field Xn of the theory El 

A = 27r( (- ^ + A + a) , Dn = 27r(Xn, 

where ( is anisotropy parameter. The principle of unitar- 
ity of the underlying field theory restricts, through the 
Kac formula, the possible values of c, A and A. For the 
2D Ising model, we have c = 1/2 and the only possible 
values are A, A = 0,1/16,1/2. Therefore there are six 


different boundary universal classes with (A, A) = (0,0) 
for periodic EC’s; (A, A) = ( 3 ^, 3 ^) for antiperiodic 
EC’s; and EC’s with (A, A) = (|4) ; ( 3 ^, O) ; (^, O) 
and ( 3 ^, ^)- Past efforts have focused mainly on peri¬ 
odic and antiperiodic EC’s. In this paper we consider 
(A, A) = ( 3 ^, 0 ), which for the Ising quantum chain is 
called duality twisted EC’s and has been considered in 
El- In what follows we will give the first exact solution 
of the 2 D Ising model with twisted EC’s at the criti¬ 
cal point. New universal amplitude ratios for finite-size 
corrections of the two-dimensional Ising model with peri¬ 
odic, antiperiodic, free, fixed and mixed boundary condi¬ 
tions have been recently presented [16]. In this letter we 
present universal amplitude ratios for the duality twisted 
EC’s universality class. 

The Ising model on a lattice with periodic EC’s and 
with a specific defect line (seam) was formulated in 
[13 di]. For the Ising model the seams are labelled by the 
Kac labels (r, 5). There are six possible partition func¬ 
tions for the Ising model with seams labelled by 

(r, s) = (1,1), (1, 2 ), (1,3), ( 2 ,1), ( 2 , 2 ), ( 2 ,3) and for the 
three of them, namely, for the seams (r, s) = ( 1 , 1 ), ( 1 , 2 ) 
and (1,3), the partition functions Z(^ s)(g) are obtained 
numerically to very high precision in [IT]- In particular 
for Z(^i^ 2 ){q) they obtain 




Xo + X1 


* I 

X_L + 


Xo +X1 


Xi, (1) 


where xo = Xo{q), Xi = X^^) and = X^iq) are 
the chiral characters and q is the modular parameter. 
The (r, s) = (1,1) and (1,3) seams reproduce the well- 
known partition function of the Ising model with periodic 
and antiperiodic boundary conditions respectively. In 
what follows we will show that the EC given by the seam 
(r, s) = ( 1 , 2 ) corresponds to duality twisted EC. 

Let us consider a square lattice rotated by 45 degrees 
(Fig. 1), in which each row has L faces and each column 
has M faces. We insert one defect line (seam) along the 
final column. The lattice thus consists of 21/ — 1 regu¬ 
lar (zigzagging) columnar edges and one defect line and 
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FIG. 1: The lattice is constructed with L faces in each row 
and M faces in each column. The line defect is inserted along 
the rightmost columnar edge (dashed line) before closing the 
lattice into a torus with periodic EC’s. 


2M regular (zigzagging) row edges. Periodic EC’s are 
imposed in both directions and in Fig. 1 this is repre¬ 
sented as identifying the light (dark) nodes on the first 
row (column) with the respective light (dark) nodes on 
the last row (column). The physical dimensions of the 
lattice, Lx and are given by 

= V2 ^ ; Ly = V2M. (2) 

The finite-size partition function for the Ising 
model can be written as = 

exp (J Y,<ij> SiSj + K Y,<ij> ^ where the first 

sum within the parenthesis is over NW-SE edges, and 
the second sum over NE-SW edges and spin variable Si 
can take the two values ±1. Since we restrict ourselves to 
the critical Ising model, we have sinh(2J) sinh(2iF) = 1. 
This condition can be conveniently parameterized 
by introducing a so-called spectral parameter u, so 
that sinh(2J) = cot(2i4), sinh(2iF) = tan(2i^), with 
0 < u < 7r/4. The anisotropy parameter is C related 
to the spectral parameter u through = sin4i4. Now 
partition function can be rewritten in the following form 

Zl.lM = Tr [T{u)f = 


where the sum is over all eigenvalues of a transfer matrix 
T{u), written as and Sn{u) = ^En(rt). 

Conformal invariance predicts that the leading finite- 
size corrections to the energies Sn take the form [ 6 ] 
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where and A^ are the conformal weights, kn^kn G N 
are label descendent levels and g is Coxeter number. 


FIG. 2: As Dynkin diagram 


The Boltzmann weights of the model are prescribed to 
the faces of a regular square lattice m 
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U I = Sl(-M)< 5 ac + Soiu) (4) 


Here a, 6 , c, d are the spin states, u {0 < u < X) is the 
spectral parameter, Sk{u) = , A = ^ is the 

crossing parameter and ^pa are the entries of Perron- 
Frobenius eigenvector of the adjacency matrix G. The 
Ising model is related to the Dynkin diagram (see 
Fig. 2 ) whose Coxeter number ^ = 4. From A^ Dynkin 
diagram we get the following allowed face configurations: 




(a,b)e{(1,1),(3,3),(1,3),(3.1)} 


We can calculate the weights of this faces by insert¬ 
ing the respective parameters in equation Obviously 
one can construct two kinds of transfer matrices U and 
V (see Fig. 3) where the rows consist of 21/ — 1 reg- 



FIG. 3: U and V transfer matrices. The rows of the transfer 
matrices consist of 2L — 1 regular faces followed by single 
seams. Single seams represents by dashed lines. 

ular faces followed by single seams. Our BC are given 
with seam (r, s) = (1,2). The indices jk take values 
1 or 3. In terms of these quantities we can construct 
the double row transfer matrix U{u)V{v) (see Fig. 4). 
Where iaGa^ ka take values 1 or 3 (a = 1 ,..., L), and over 
the variables /ci, /c 2 ,..., a summation is performed. On 
other vertexes assume non-fluctuating heights equal to 
2 . The double row transfer matrix U{u)V{v) satisfies 
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FIG. 4: The UV double row transfer matrix consists of L — 1 
blocks and the last dashed block. 


the functional equation 

U{u)V = cos I -\-i ( — 1)^ sin R 

where / is the identity matrix, and R is a square matrix 
with anti-diagonal entries equal to 1 with remaining en¬ 
tries 0 [both are (2^ x 2^) matrices]. For the eigenvalues 
A{u) of the double row transfer matrices UV of the Ising 
model with duality twisted EC’s we can obtain 

L-l 

A = e2“ Y[ (e2» sin cos ^k )' (5) 

k=l 

where u = u — cpk = 2 ( 2 L-i) ~ 

chosen arbitrarily so we have 2^~^ eigenvalues. The re¬ 
maining eigenvalues can be found by taking the complex 
conjugate of (§. From Eq. (§ one can easily obtain 
the expressions for the absolute value of the eigenvalue 
|A| and the argument 0. Let us now consider the largest 
eigenvalue Aq, which corresponds to the case when all 
/i/c = 1. The derivation of the asymptotic expansion of 
logAo = log|Ao| + iOo can be divided into two parts. 
First, with the help of the Euler-Maclaurin summation 
formula we can derive the asymptotic expansion of the 
logarithm of the absolute value of Aq 
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log |Ao| = Lxfoo ^ ^ 

/c=0 


2 ^( 2 A: + 2 )! 



where ^ 2 / 0+2 are the Bernoulli numbers, f{x) = 
^ log [1 + sin(4i4) sinx] and foo = —^fQf(x)dx. Next, 
with the help of the Boole summation formula [18] the 
asymptotic expansion of the argument Oq can be written 
in the form 


^2n+l (0) / TT \ 

2-(2n + l)! 


where g(x) = i arctan ^ 2 n+i( 0 ) are the 

Euler polynomials. Thus the leading finite-size cor¬ 
rections to the ground state energies £0 = = 


— ^ log Aq take the form 



On the other hand, £0 given by Eq. @ Comparing Eq. 
^ with Eq. © for the Ising model with c = 1/2 and 
^ = 4 we see that A = ^ and A = 0. Thus we have 
shown that duality twisted EC’s with (A, A) = ( 1^5 O) 
correspond to the BC with seam (r, s) = (1,2). 

Let us now consider the other eigenvalues given by Eq. 
(§, which correspond to some combination of the /i/^. 
Denote by Ajpj. the eigenvalue with fip = —1 and the re¬ 
maining Ilk having the value /i/c = +1. The leading finite- 
size corrections to the exited state energies £p = = 

-^logA{p} and £l-p = ^^L-p = -^logA{i_p} 
can be written in the form 


£p 


£l 


-p 
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where is given by Eq. (©• Comparing Eqs. §, ® 
with Eq. (§ for the Ising model with c = 1/2 and ^ = 4 
we can see that for the exited state £p we have A = 1/16, 
/cn = 0 and A = l/2, 1 and for the exited state 

£l-p we have A = 1/16, = P and A = 0, = 0- 

Now we have all the necessary information to start the 
calculations of partition function for the Ising model with 
duality twisted EC’s Eor large and Ly (always 

keeping the ratio LyjLx constant) we have 


Z Xjy 




where ^( 1 , 2 ) (^) is the universal conformal partition func¬ 
tion. Let us now find the general form of the eigenvalues 
with significant input in the partition function. Erom 
(§ we can see that these are eigenvalues for which al¬ 
most all /i/c = 1 and some jik are allowed to take the 
value —1 only \i k L oi L — k L. Any ’’signifi¬ 
cant” eigenvalue will be specified by two sets of indexes 
K = {ki, fe, • • • , and K = {^ 1 , fe, * * * , where 
ki = L — ki with ki < k 2 < ks <—'< km ^ L and ki < 
< ^3 < • • • < ^ A so that Ilk. = — 1 ,/i/ci = 

the other p’s are +1. Thus from above it is easy to get 

Z(i,2){q) = + ( 9 ) 

KK 


_ ^TviLy —4iu 

where q = e is a modular parameter. Eor 

further calculations it is convenient to introduce the oc¬ 
cupation numbers Sk = and Sk = • la terms of 

these quantities the universal conformal partition func- 























tion can be rewritten as 




^ ^ q 2°4+2)^fe 


= q 2°4+16gr 2°4 (^1 + grfe) jj 2^. 

k=l k=l 


Taking into account the expressions for the chiral charac¬ 
ters Xo, Xi, HH and adding the conjugate part of the 
eigenvalue set, we can get the universal conformal parti¬ 
tion function given by Eq. §■ Thus we have obtained 
analytically the universal conformal partition function 
for the Ising model with duality twisted EC’s ^( 1 , 2 ) (^) 
which confirms the numerical result of m- 

Now we will present the new set of the universal 
amplitude ratios. Let us denote the free energy per 
spin /, the inverse correlation lengths and 


of our critical Ising model as = — ^log(|Ao|), 

(i' = t1»s(| 4^|) “d (li, = 

We find that subdominant finite-size corrections to scal¬ 
ing should be to the form aklL‘l^~^ for the free energy / 
and bk{p)I and Ck{p )/for inverse correlation 


lengths ^ and respectively, where the coefficients 

bk{p) and Ck{p) can be written in the following form 


dk 


hip) 

Ck{p) 


2 *^-i( 2 fc)! 

^ 2 fe-l( 2 p_ l) 2 fe-l 

2 fe- 2 ( 2 )!c- 1 )! 
1 ^ 2 fc-l„ 2 fe-l 


2k 

j( 2 fe-l)^Q) 


(2fc-l)! 




The coefficients ai, 6 i(p) and ci(p) are universal and re¬ 
lated to the conformal anomaly number (c), the confor¬ 
mal weights (A, A), and the scaling dimensions of the p- 
th scaling fields of the theory. The coefficients a/e, bk{p) 
and Ck{p) for k > 2 are non-universal, but ratios of these 
coefficients rp{k) = and rL-p{k) = are uni¬ 

versal and given by 


Tpik) — 


4.k{2p-lf^-^ 




2k 


rL-p{k) = 


Ak{2p) 


2k-l 




2k 


Thus we have obtained a new set of the universal am¬ 
plitude ratios rp{k) and rL-p{k). The case for k = 1 
is trivial, since rp(l), rL-p(l) are the ratios of universal 
coefficients ai, bi{p) and ci(p). The case k > 2 is non¬ 
trivial. Below we will show that the ratios rp{2) and 
rL-p{2) which are given by 


rp{2) = -240(2p-l)^ (10) 

rL-pi2) = -1920p3 (11) 


can be obtain from conformal field theory. The finite- 
size corrections to Eq. © can in principle be com¬ 
puted in perturbative conformal field theory. In gen¬ 
eral, any critical lattice Hamiltonian will contain correc¬ 
tion terms to the fixed-point Hamiltonian H = (He + 


^k9k /_£ /2 where pk is a non-universal con¬ 

stant and (j)k{v) is a perturbative conformal field wit^ 
scaling dimension Xk- To the first order in the pertur¬ 
bation, the energy gaps {£n — ^o) and the ground-state 
energy {£q) can be written as 


£n — £0 


£0 


27r 


C^n H~ 27r ^ ^ pk {£^nkn ^O/co) 




fo,c + 27r E 9kCokO 

k 




where Cnkn are universal structure constants. In the case 
of the cylinder the spectra of the Hamiltonian are built 
by the irreducible representation A, A of two commut¬ 
ing Virasoro algebras Ln and Zn- The leading finite- 
size corrections (l/L^) can be described by the Hamilto¬ 
nian with a single perturbative conformal field (j)i{v) = 
L‘^_ 2 {v)^L‘^_ 2 {v) with scaling dimension = 4 [19]. Thus 
the ratio rn( 2 ) are indeed universal and given by 


rn( 2 ) = 


Cnln — Ci 


010 


C'olO 


The universal structure constants Cnin can be obtained 
from the matrix elements (n|(/)i( 0 )|n) = (2'kILxY^ Cnin 
[ 20 ] , which for non-degenerated states have already been 
computed by Reinicke m- 


Cnin — 


(A + r) 

(s)“- 


12 (A + 1)(2A + 1) ) ^ ’ 

^ ^ [^"(5c - 8 ) - (5c + 28)] 5a, 0 - 


Let us consider the case c = 1 / 2 . Eor the two-dimensional 
Ising model with duality twisted EC’s the ground state 
|0 >, the excited state \p > and the excited state \L—p> 
are given by |0 >= |Ao = r = 0 ; Aq = 0 , f = 0 >, 
\P >= |Ap = ^,r = 0; Ap = = p-l > and \L-p >= 

\^L-p = ^,r=p; Al-p = 0 , f = 0 >. 

Eor non-degenerated states the universal structure con¬ 
stants Cnin (n = 0,p and L — p) can be obtained from 
Eq. (12) and for ratios rp{2) and rL-p{2) one can obtain 


rp{2) = -240(2p-l)3 
rL-p( 2 ) = -^^pilQp^ 


(13) 

(14) 


■3p-2). 


Thus we can see that Eq. ( |13[ ) coincides with Eq. ( |1Q[ ) 
for all values of p, while Eq. (14) coincides with Eq. 


(11) only for p = 1 and 2. The excited states jL — p > for 
p > 3 are degenerated and one cannot apply the Reinicke 
formula (12). Eor degenerated states the calculations of 


the universal structure constants Cnkn is not straight¬ 
forward, but for the case of the Ising model it can be 
done. Eor the Ising model we have calculated the univer¬ 
sal structure constants Cnkn for degenerated states and 
find that the results are in complete agreement with Eq. 
© for all values of p. 
























In this Letter we considered the Ising model with du¬ 
ality twisted EC’s. We have showen that EC’s with 
(A, A) = (^,0) correspond to the seam (r, s) = (1,2). 
We derive exact expressions for all eigenvalues of the 
transfer matrix for the critical Ising model with the du¬ 
ality twisted EC’s. We derive explicitly the universal 
conformal partition function ^( 1 , 2 ) (^) and show that it is 
consistent with numerical results of [T7] . We find that the 
ratio of the subdominant finite-size corrections to scaling 
in the asymptotic expansion of the free energy / and 


the inverse correlation lengths and are univer¬ 

sal and give their exact value. We verify this universal 
behavior using a perturbative conformal approach. ^ 
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